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ABSTRACT 
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CNJ , Deep spectroscopy of the two millijansky radio galaxies LBDS 53W069 and LBDS 

■ 53W091 has previously shown them to have old (> 3 Gyr) stellar populations at 

z ~ 1.5. Here we present the results of Hubble Space Telescope (HST) observations 
with the Wide Field and Planetary Camera 2 (WFPC2) in F814W and with the 
I ■ Near-Infrared Camera and Multi-Object Spectrograph (NICMOS) in F110W. We find 

that 53W069 has a de Vaucouleurs r 1 / 4 profile in both the F814W & F110W data 
£3 . with a mean effective radius of 0'.'30 ± 0'.'06 (2.7 ± 0.5 kpc). The restframe U - B 

colour gradient is consistent with that of present-day ellipticals, requiring a stellar 
population of super-solar (3Zq) metallicity that formed on a very short timescale at 
high rcdshift (z > 5). 53W091 has a regular r 1 / 4 profile in F110W with an effective 
radius of 0'.'32 ± 0'.'08 (2.9 ± 0.7 kpc). The F814W profile is more extended and is 
5_i ■ consistent with the presence of a blue exponential disk that contributes 20 ± 10% of 

the flux within r e . We find a restframe U — B colour gradient which is significantly 
larger than that observed in field ellipticals at z ^ 1, implying a stellar population of 
mixed metallicity (1-3Zq) that formed in a high-redshift rapid burst. 

We have compared these two LBDS radio galaxies with the Kormendy relations 
of ten 3CR radio galaxies at z ~ 0.8 and a sample of cluster ellipticals at z ~ 0.4. 
The LBDS galaxies follow the Kormendy relation for the more radio-luminous 3CR 
galaxies, assuming passive evolution of their stellar populations, although they are 
smaller than the 3CR galaxies whose mean effective radius is 12 kpc. Their sizes and 
radio luminosities are consistent with scaling relations applied to the 3CR galaxies, 
in which both radio power and effective radius scale with galaxy mass. Compared 
with the sample of cluster ellipticals, 53W069 & 53W091 lie well within the scatter of 
the Kormendy relation. We conclude that the hosts of these millijansky radio sources 
at z ~ 1.5 are passively-evolving elliptical galaxies that will evolve into ordinary L* 
ellipticals by the present day. 

Key words: galaxies: active — galaxies: elliptical and lenticular, cD — galaxies: 
evolution — galaxies: individual (LBDS 53W069, LBDS 53W091) 
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1 INTRODUCTION 



factor of ~ 200 fainter than the 6C survey, t hus probing 



How did elliptical galaxies lorm and evolve? ihere are two 



lower radio luminositie s and higher redshifts. Windhorst 



mecharilymy iliui have been proponed, first is Hie moiroritttc 

collapse model, in which elliptical galaxies f orm at high red 



shifts in an intense bu rst of star formation (Eggen, Lynden 



Bell & Sandage 1962). Essentially all the final mass of the 
galaxy is already present within its potential well at the time 
of formation. The second mechanism is the merger model, in 
which ellipticals form at relatively lo wer redshifts from the 
merg er of many smaller galaxies (e.g. 



Kauffmann & Chariot 



Keel & Pascarelle (1995 ) showed that the radio galaxy LBDS 
53W002 has a weak AGN and a dominant r 1//4 profile with 
the colours of a ~ 3 x 10 8 year young stellar population at 
z = 2.39, suggesting that these weak radio sources could 
have formed rather mature ellipticals at high redshifts. Two 
more of these sources have proved particularly useful for 
the study of galaxy evolution, due to the deep spectra that 
we obtained with the Keck telescope. LBDS 53W091 has 
a redshift of 1.552, and its restframe UV spe ctrum is best 



199£j). In this model, the mass of the galaxy increases with mode l ed by a stellar popul ation > 3.5 Gyr old ( |Dunlop et al 



time asl more and more smaller galaxies are cambalized by 
the forrnihg elliptical. 



199C; gpinrad et al. 1997). LBDS 53W069 is at z = 1.432 



and has a best-fitting stellar population of age ^ 4 Gyr ( Dey 



Re ;ent observations suggest that both mechanisms may 199 7 ■ D un lop 1999b). Spinrad et al. (1997) investigated the 



have a role to play, lor example, pJllis et al. (199Y ) lound 
that spheroidal galaxies in clusters at z ~ 0.5 required a for- 
mation redshift (by which we mean the redshift at which the 
dominant s tellar population formed) of zj ^ 3. In contrast, 
Zepf (1997) found that the small number of very red galaxies 
in deep optical/infrared surveys suggested either that ellip- 
ticals form at moderate redshifts (and are enshrouded by 
dust) , or that they assemble through the merging of smaller 
galaxies. Similarly Pascarelle et al. (1996) identified a group 
of Lyman-a emitters at z — 2.39 that they suggested may 
subsequently merg e at a lower redshift i nto one or more 
luminous galaxies. Treu & Stiavelli (1999) concluded from 
their study of NICMOS parallel fields that a significant frac- 
tion (10-66%) of the elliptical galaxy population formed at 
Zj > 3, but the rest may have been formed (or at least 
assembled) at lower redshifts. 

One important method of studying high-redshift (z > 
1) ellipticals is via radio selection. At low redshifts, lumi- 
nous radio sources are almost exclusively hosted by giant 
elliptical galaxies containing old stellar population s (e.g 



Kron, Koo fc Windhorst 1985; Taylor et al. 1996 



Nolan 



et al. 2001a). To the extent that radio galaxies at high red- 
shifts are also hosted by giant ellipticals, they can be used 
to study the evolution of the elliptical galaxy population. 
Lilly fc Longair (1984 ) obtained infrared K magnitudes for 
a subsample of the 3CR radio survey, from which they con- 
structed a Tf-band Hubble diagram. They concluded from 
this K—z relation that luminous radio galaxies at z ~ 1 are 
giant ellipticals w ith passively-evolving stellar p opulations. 
Lilly et al. (1985| ) and |Dunlop fc Peacock (199o[ ) found that 
the K—z relation for less-powerful radio galaxies led to a 
similar conclusion. 



At redshifts z > 0.6, Eales et al. (1997) found that the 
6C radio galaxies, with radio luminosities a factor of ~ 5 
lower than the 3CR, were also on average 0.6 mag fainter 
than 3CR in the if -band. Roche et al. (1998 ) investigated 
the A"'-band morphologies of ten of these 6C sources at z ~ 
1.1, showing that seven of them were normal ellipticals and 
the other three were ongoing or recent mergers. The radii of 
the 6C galaxies were significantl y smaller than those of the 



3CR sources at similar redshift s (Best et al. 199£ ; Zirm et al 



199£; McLure & Dunlop 2000), indicating that their fainter 



K magnitudes were due to their smaller size, not solely due 
to the difference in power of their AGN. 

The Le i den-Berkeley Deep Surve y (LBDS; Windhorst 
ct al. 1984; Waddington et al. 200C), with a flux density 



limit of 1 mjy at 1.4 GHz, is sensitive to radio sources a 



age determination in detail and showed how the minimum 
age of the universe at z — 1.5, as required by the age of 
53W091, placed limits on the all owed values of the co smo- 
logical parameters (Ho, Qm, Oa)- Stockton et al. (1995) sim- 
ilarly showed that the age inferred from the Keck spectrum 
of the radio galaxy 3C65 (~ 4 Gyr at z = 1.175) required a 
formation redshift of Zf ^ 5. 

Several authors have qu estioned the reliability of these 
age estimates, for example Bruzual & Magris (1999) and 
Yi et al. (200C|) deriv e ages of 1.5-2 Gyr for 53W091. How- 
ever, punlop (1999^ ) argued that such young ages are only 
deduced if the near-infrared photometry is included in the 
model fitting; if the fitting is confined to the spectroscopic 
data then a variety of stellar population sythesis codes con- 
sistently produce ages of ^ 2.5 Gyr (see also Nolan et al. 
2001b). We will not discuss this age controversy further here, 
but simply note that nothing in the current paper depends 
crucially on the precise ages of these sources. 

In this paper we will investigate the morphologies of 
these two radio galaxies with the Hubble Space Telescope's 
(HST) Wide Field and Planetary Camera 2 (WFPC2) 
and Near-Infrared Camera and Multi-Object Spectrograph 
(NICMOS). At z ~ 1.5, the 4000-A break is straddled by the 
F814W and F110W filters, and thus the emission observed 
through these two filters is dominated by the young and old 
stellar populations respectively. In section 2, we describe the 
observations and discuss the processing steps that were ap- 
plied to the data. In section 3, we investigate the surface 
brightness profiles and colour gradients of the two sources. 
The location of the galaxies on the Kormendy relation is the 
topic of section 4 and finally, in section 5, we comment on 
the apparent relations between radio luminosity, black-hole 
mass and effective radius of these sources. Our results are 
summarized in section 6. 

We use AB magnitudes unless otherwise noted, and 
denote magnitudes in the three HST filters by 7si4 for 
F814W (WFPC2), Jno for F110W (NICMOS) and H 160 
for F160W (NICMOS). We assume a flat cosmology with 
H = 65 km s" 1 Mpc" 1 , Q M = 0.3, Ov = 0.7. 



2 OBSERVATIONS AND DATA REDUCTION 

In 1997 August and 1997 October, we observed 53W069 and 
53W091 with WFPC2 and NICMOS on the HST. 53W069 
was observed for 2 orbits (a total integration of 5100 s) in 
F814W with WFPC2, and for 3 orbits (7680 s) in F110W 
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Figure 1. WFPC2 F814W and NICMOS F110W images of 53W069. The logarithmic grey-scale has been chosen 
so that both images have the same range in flux density (Fx). Each image is 3 X 3 arcsec 2 . The F814W image has 
been rotated in order to display it at the same orientation as the F110W image. 



with NICMOS camera 2. 53W091 was observed for 2 orbits 
(5300 s) in F814W with WFPC2, and for 2 orbits (5120 s) 
in F110W with NICMOS camera 2. Each orbit devoted to 
WFPC2 imaging was split into two exposures, and consec- 
utive orbits were offset by 0'.'35. In each orbit devoted to 
NICMOS imaging, five 512 s exposures were taken using a 
spiral dither pattern with l'.'91 offsets. In 1998 August, we 
observed 53W091 and its companions for 10 orbits (25 600 s) 
in F160W with NICMOS camera 2. We independently anal- 
yse tho se data and describ e the results in a complementary 
paper (Bunker et al. 2002). 



Table 1. Photometry 



The WFPC2 exposur es were re duced using the stan- 
dard pipeline processing ( Voit 1997 ). For each target, the 
four dithered exposures were com bined into a single im- 



age using the drizzl ing technique (Hook & Fruchter 1997 



Fruchter et al. 1997 ). A pixel size of 0'.'05, one-half of the 
original pixel scale, was used, with a drop size of 0.6 (the 
parameter PIXFRACT). Offsets were calculated from the 
cross-correlation of ~ 10 bright sources in each exposure; 
no rotation was required. The drizzled images are shown in 
Figs | & | 

We reduced the NICMOS imag es using a combi nation of 
processing techniques from STSDAS (Dickinson 1999) and our 
own methods. The raw images were partially processed with 
the standard pipeline task CALNICA, using the best refer- 
ence files available from the STScI. We tried using the latest 
temperature-dependent darks in order to remove the shad- 
ing, however we obtained better results with the standard 
dark reference files and used these instead. Non-linearities 
in the bias level of each image quadrant were removed with 
BIASEQ. Bias jumps were not corrected for in the data, as 
only one image was significantly affected and BIASEQ was not 
successful in removing the jumps. Cosmic rays were then re- 
moved by continued processing with CALNICA. The variable 
bias level from quadrant to quadrant, or "pedestal" , was cor- 
rected using the pedsky task, which also subtracted the sky 



Source 


Total AB 


magnitude 


Colour* 




^814 


JllO 


^814 — JllO 


53W069 


24.19 ±0.03 


22.45 ±0.02 


1.72 ±0.05 


53W091 


24.14 ±0.03 


22.56 ±0.03 


1.57 ±0.05 



*The colours are defined in a circular aperture of 2" diameter. 



background from the images. Then the data were flattened 
using CALNICA. 

We applied a further correction to the data, adapted 
from standard ground-based infrared imaging methods. 
Since the NICMOS camera consists of four physically sep- 
arate sub-arrays, we divided each partially-reduced image 
into its four separate quadrants. For each target, all the ex- 
posures of each quadrant were stacked and a median image 
calculated. Given that the exposures were dithered, this pro- 
duced a map of the residual instrumental features (devoid 
of any astronomical objects) that were not removed by the 
flat-field reference file. These "super-sky" images for each 
quadrant were scaled to the mean of all four quadrants in 
order to preserve the quadrant-to-quadrant photometric ac- 
curacy. We then tried correcting each exposure by either (i) 
dividing by this super-sky after scaling it to have a mean 
of unity, or (ii) subtracting the super-sky from each expo- 
sure. For each object, we chose the method that produced 
the best image (i.e. the lowest rms noise). For the F110W 
images of both targets, subtracting the super-sky was most 
successful, suggesting that the instrumental features left af- 
ter CALNICA processing were most likely due to incomplete 
dark subtraction. In particular, this removed the residual 
amplifier glow from the corners, while making a negligible 
difference to the center of the image, where the target was 
located (thus photometric accuracy had not been compro- 
mised by residual dark current beneath the target) . Finally, 
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Figure 2. WFPC2 F814W and NICMOS F110W images of 53W091. The logarithmic grey-scale has been chosen 
so that both images have the same range in flux density (Fx). Each image is 3 X 3 arcsec 2 . The F814W image has 
been rotated in order to display it at the same orientation as the F110W image, resulting in the apparent blurring 
of the drizzled data. 



all the exposures were combined, on a quadrant by quad- 
rant basis, to produce a mosaic for each object (Figs |l| & |^). 
The data were calibrated using the most recent photometric 
calibration parameters from the STScI. 



Th e imag e detection program SEXTRACTOR (Bertin fa 



1996) was used to calculate both aperture and to- 



Arnout 

tal ( "mag_auto" ) AB magnitudes for the two radio galaxies 
in each of the filters. The results are given in Table [j]. Al- 
though the total magnitudes give the best measurement of 
the flux detected in each individual image, the colours of the 
sources are properly defined in fixed apertures, thus we give 
the colours in Table ^ for 2" diameter circular apertures. 
We note that integration of the best-fitting model surface 
brightness profiles (see section 3.1 below) gives magnitudes 
consistent with these measurements. 



3 MORPHOLOGIES 

3.1 Surface brightness profiles 

Surface brightness profiles were extracted for the two radio 
galaxies by fitting elliptical isophotes to the sky-subtracted 
images (Figs | fa |ij) with iraf's ellipse routine. We fixed 
the centre of each isophote and allowed the position angle 
and axial ratio to vary freely. Model surface brightness pro- 
files were extracted in the same way from two-dimensional 
exponential and de Vaucouleurs r 1//4 models convolved with 
a model point -spread function (PSF) from Tiny Tim ( Krist 
fa Hook 1999). The models were subsampled by a factor of 
1, 2 or 5, convolved with a subsampled PSF and then resam- 
pled to the (drizzled) WFPC2 and NICMOS pixel size, be- 
fore being extracted to a one- dimensional profile. The mod- 
els were fitted to the data for r > 1 pixel, corresponding 
to 0'.'05 for WFPC2 (drizzled) and 0'.'075 for NICMOS - at 
smaller radii, the fit was dependent on the degree of sub- 
sampling of the PSF. 
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Figure 3. Surface brightness profiles along the major axis 
(points) and best-fitting models (lines) for 53W069. Both profiles 
follow a de Vaucouleurs i- 1 / 4 law, with effective radii of r e = 0'.'33 
in F110W and r e = 0'.'28 in F814W. 



The best-fitting models for 53W069 are de Vaucouleurs 
r 1/4 profiles in both the F814W fa F110W filters (Fig. |). In 
F814W we find an effective radius of r e = 0"28 ± f/.'08, and 
in F110W the profile is very similar with r e = C/.'33 ± tf/09. 
The reduced chi-square for these models is \r — 2, com- 
pared with a value of — 10 for the best-fitting exponen- 
tial models. The weighted mean of the F814W fa F110W 
effective radii is r e = 0'.'30 ± 0'.'06 or 2.7 ± 0.5 kpc. Morpho- 
logically, 53W069 is thus a relaxed elliptical galaxy, fully 
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Figure 4. Surface brightness profiles along the major axis 
(points) and best-fitting models (lines) for 53W091. The F110W 
profile follows a de Vaucouleurs r 1 / 4 law, with effective radius of 
r e = rX'32. The F814W profile is the sum of an r 1 / 4 profile with 
effective radius r e = 0732 plus an exponential with scale-length 
ro = 075 contributing 20% of the flux within r e . 



consistent with its spectroscopic identification as an old, 
passively-evolving elliptical. 

For 53W091, the best-fitting model in F110W is a de 
Vaucouleurs r 1/4 profile with r e = C/.'32±0'.'08 (Fig. §). This 
model has a reduced chi-square value of \ 2 r — 1-8 compared 
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Figure 5. / 814 - J 110 colour gradients for 53W069 & 53W091. 
Crosses denote data within the central two pixels for which pixe- 
lation effects and differences in the WFPC2 and NICMOS PSFs 
contribute as much as 0.1-0.3 mag to the colour gradient. The 
horizontal lines are the colours predicted by a high-redshift burst 
of star formation, for two different metallicities: solar (dot-dash 
line) and 3Zq (dashed line). 



with a value of £ = 7.3 for the best exponential model. In colours of the galaxy, w h ich become bluer with increasing 



Bunker 



ot al. (2002) wo chow that tho F160W profilo eim 



ilarly fallows a Hp Vaiimnlpurs law with r a — 0"3? + 0"03 

For the F814W profile, the best-fitting single component 
model is an r 1 / 4 profile with r e — 0'.'6 ± 0'.'2, for which 
X 2 , = 3.3 (compared with \ 2 r = 18 for the best exponen- 
tial model). Given the large difference between the effective 
radius of this best-fitting model and that of the F110W & 
F160W models, we also fitted the F814W profile with a two- 
component model. We used the best-fitting F110W/F160W 
de Vaucouleurs model with r e = 0'.'32 and added an expo- 
nential disk, varying both the scale-length of the exponential 
and the relative flux of the two components. The best-fitting 
composite model (with xt ~ 2.6) has a disk with scale- 
length r = 0'.'5 ± 0'.'2 (4 ±2 kpc), that contributes 20±10% 
of the flux within the half-light radius (Fig. At z — 1.552, 
the F814W filter measures the galaxy's restframe emission 
below the 4000- A break, and is thus sensitive to bluer stars 
than F110W. Such a population will make an even less sig- 
nificant contribution to the flux at longer wavelengths, thus 
it is quite consistent that we do not see a disk component 
in the F110W data. 



3.2 Colour gradients 

The stellar population of an elliptical galaxy is not uniform, 
but varies with radius. This is seen as a gradient in the 



radiu s (Franx et al. 1989 



Peletier et al. 199C ; Famura et al 



200C| ). If 53W069 and 53W091 are the passively-evolving 
precursors to present-day ellipticals, we should expect to see 



Peletier et al. (199C) 



such characteristic colour gradients 
show that the U — B colour gradient of nearby ellipticals 
is —0.1 mag arcsec -2 per dex in radius. The I%n and Jno 
bands of our radio galaxies at z ~ 1.5 approximately corre- 
spond to rest-frame U and B, t hus we should expec t to see 



Peletier et al. (199(]| ). 



a gradient of the same order as 

In Fig. ^, we plot the colours of the two radio galaxies as 
a function of radius. Small differences between the (drizzled) 
WFPC2 and NICMOS PSFs dominate the colour gradients 
within the central 1-2 pixels. Specifically, the radii contain- 
ing (50%, 90%) of the encircled energy are (0'.'08, 0'.'37) and 
(0708 , 0748) for the WFPC2 F8UW and NICMOS F110W 
PSF models respectively. Examination of our PSF-convolved 
galaxy models suggests that these differences result in arti- 
ficially redder colours within the central pixels, contributing 
as much as 0.1-0.3 mag to the colour gradient at r < 0715 
(crosses in Fig. |E|). At larger radii this effect is less than 
0.1 mag, so we restrict our discussion to r > 0715. Given 
the limited quality of the data and the absence of PSF stars 
in the NICMOS images, we decided there was little to be 
gained from attempting detailed PSF matching to better 
than 10%. 

It can be seen in Fig. |^ that the errors in the 7si4 — 
Jno colours of 53W069 are too large to measure a colour 
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change of —0.1 mag arcsec" 2 across the extent of the galaxy 
(less than 1 dex in radius). Thus the colours of 53W069 
do not preclude a gradient comparable to that of present- 
day ellipticals, but they are also consistent with there being 
no colour gradient. Windhorst et al. (1994, 1998) similarly 
found little or no colour gradients in several other LBDS 
radio galaxies and field galaxies at 0.3 < z < 2.4 observed 
with HST. In particular, any colour gradient in the elliptical 
radio galaxy 53W002 at z = 2.39 must be smaller than 



outer regions (Hinkley & Im 2001), although the observed 



—0. 3 m ag across the 1" (9 kpc) image ( Windhorst et al, 
199 8j). fo ie (restframe) U — B gradients of early-type cluster 
and field galaxies at z ~ 0.4-1 have been found to typically 
span —0.2 to 0.1 m ag per dex in radius ( Baglia et al. 2000; 
Hinkley fc Im 2001), again consistent with 53W069. 



In contrast, there is a significant gradient in the 7gi4 — 
Jno colour of 53W091. This gradient is in the same sense as 
that seen in nearby ellipticals (i.e., the galaxy is redder in 
the center and becomes more blue at larger radii), but it is 
much greater in magnitude. The best fit to the colours gives 
a gradient of — 1.0±0.2 mag arcsec -2 per dex in radius. This 
clearly indicates the presence of an extended blue component 
in the galaxy's morphology, and is consistent with the results 
of our model fitting above. There appears to be a "jump" in 
the colour gradient at r ~ Of! 4 (see also the F110W profile), 
but this is not significant given the errors, and the data 
are consistent with a constant gradient out to r ~ Cf.'6. We 
note that the position angle (PA) of the major axis of the 
optical/infrared emission (~ 45 degrees) is perpendicular 



to t he P A of the radio lobes (131 degrees; Spinrad et al, 
199? |), | o it is unlikely that this extended blue emission is 
due to the radio-optical alignment effect. We emphasize that 
the difference in WFPC2 and NICMOS PSFs is unlikely to 
have produced this gradient because: (i) it is not seen in the 
data for 53W069; and (ii) the profile fitting above includes 
a convolution with the PSFs, and clearly shows that the 
galaxy is more extended (i.e. has a larger scale-length) in 
the F814W data than in the F110W data. 

We compared the observed colour gradients of 53W069 
and 53W091 with the predictions of single collapse multi- 
zone models of e lliptical galaxies (Menanteau, Jimenez & 
Matteu:ci 20011). We used an IMF with lower and upper 
cutoffs of 0.1 Mq and 120 Mq respectively, and a range of 
formation redshifts. None of the models could reproduce the 
observed gradients, demonstrating that the very red cores 
found in these two galaxies are inconsistent with an extended 
period of star formation in their centres. We then modified 
the models to force the star formation to stop after 0.05 Gyr. 
The resulting colours for two different metallicities, solar 
(Zq) and 3Zq are shown in Fig. |f| We find that the models 
fit the observed gradients if star formation took place on a 
short time scale in both galaxies and if they were formed 
at redshift z > 5. 53W069 is best fit by a population that 
has, on average, supersolar metallicity (3Zq). 53W091 has 
a core that is metal-rich (3Zq) and at larger radii (r > 0'.'4) 
is better fit by a solar metallicity population. 

These results are in excellent agreement with those 
of |Nolan et al. (2002j) who find that the UV spectra of 



these galaxies are similarly best-fit with populations of solar 
(53W091) and supersolar (2Zq, 53W069) average metallic- 
ity. Studies of z ~ 0.4-1 ellipticals similarly indicate that the 
colour gradients of early-type galaxies are due to metallicity 
gradients, falling from ~ 2.5Zq at the centre to solar in the 



gradients are not as steep as that in 53W091. 

As discussed above, the best fit to the F814W profile of 
53W091 is obtained by adding a disk-like exponential com- 
ponent to the best-fitting de Vauc ouleurs profile m easured 
in the F110W and F160W bands. |Dunlop (1999a| ) showed 
that the difference between the restframe UV spectrum of 
53W091 and that of 53W069 can be fitted by a blue stellar 
(F0V) spectrum, contributing 10-20% of the UV flux. Sim- 
ilarly, we have shown that ~ 20% of the restframe [/-band 
flux within the half-light radius (r e = 0'.'32) could originate 
in a disk component. With a 1" slit and seeing of O'.'8-l", the 
spectroscopy of 53W091 was dominated by the flux within 
r e , and we find consistency between the fraction of excess 
blue light in the spectroscopy and the disk component in the 
imaging. This extended blue component is also consistent 
with having lower metallicity than the core of the galaxy. 

Combining these results leads us to conclude that 
53W091 is either: (i) an elliptical galaxy with an old, 
metal-rich, passively-evolving stellar population (just like 
53W069), surrounded by a faint disk of younger or lower- 
metallicity stars; or possibly (ii) an elliptical galaxy with an 
old, passively-evolving stellar population whose metallicity 
decreases with radius. We note that the age of the dominant 
red stellar population in 53W091 would have been somewhat 
underestimated du e to the presence of a you nger stellar com- 
ponent in the disk. Menanteau et al. (2001a| ) find that about 
half of field ellipticals at z ~ 1 must be undergoing recent 
episodes of star formation, due to their irregular internal 
colours and blue cores. The two radio galaxies 53W069 & 
53W091, in contrast, belong to the other 50 per cent of field 
ellipticals that have evolved passively since very high red- 
shifts. 



4 THE KORMENDY RELATION 

A key characteristic of dynamically relaxed el liptical galax- 
ies is the existence of the Funda mental Plane (Djorgovski & 
Dressier et al. 1987 ) - a scaling relation between 



Davis 1987 



the luminosity L, effective radius r e and velocity dispersion 
a of ellipticals: L oc r^a b . Given the difficulty of measuring 
velocity dispersions, particularly at high redshift, a more 
useful relation is the projection of the Fundamental Plane 
onto the effective radius -surface brightne ss (fj, e ) plane, or 
the Kormendy relation ( Kormendy 1977 ). Cluster ellipti- 
cals (|Bar gcr et al. 1998 ; Zicglcr ct al. 1999[), field ellipticals 
(pchade et al. 1996 ) and radio galaxies (McLure & Dunlop 



200C ) have all been shown to follow the Kormendy relation 
out to z ~ 0.5-1.0, assuming passive evolution of their stellar 
populations. Here we consider if our two high-redshift LBDS 
radio galaxies are consistent with these results. At z ~ 1.5 
the F110W filter approximately samples the restframe B- 
band. We thus compare the F110W data on 53W069 and 
53W091 (Table |^) with two samples of elliptical galaxies 
observed with HST in restframe B. 



4.1 Comparison with 3CR radio galaxies at 
z ~ 0.8 

In Fig. ^, we plot the surface brightness at the effective ra- 
dius /x e against the effective radius r e for the two LBDS ra- 
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Table 2. Data for the F110W (restframe B-band) Kormendy 
relation. fi e is the surface brightness at the effective radius (r e ), 
and (fi e ) is the mean surface brightness within r e . 



Source 




Me 


{He) 




(arcsec) 


(mag arcsec -2 ) 


(mag arcsec -2 ) 


53W069 


0.33 ±0.08 


22.6 ±0.1 


21.7 ±0.1 


53W091 


0.32 ±0.08 


22.4 ±0.1 


21.5 ±0.1 



dio galaxies and a sample of 3CR radio galaxies at z ~ 0.8 
(iMcLure & Dunlop 2000h. The F110W and F814W filters 



both sample the rest-frame B-band at z — 1.5 and z = 0.8 
respectively. We have thus corrected the Jno-band surface 
brightness of 53W069 & 53W091 to their equivalent 7 8 i 4 - 
band values at z = 0.8, minimizing any error introduced by 
the possible wavelength dependence of r e . We made three 
corre ctions derived using the spectral sythesis models of 
both Fioc fc R occa-Volmcrari fie (1997 ) and the 1996 revi- 
sion of Bruzual & Chariot (1993). We used their elliptical 



and burst models respectively, and refer the reader to those 
papers for technical details such as IMFs. The F110W fil- 
ter is significantly wider than F814W, thus requiring a K- 
correction of ±0.6 mag to the F110W data. Second, pas- 
sive evolution of ±0.7 mag from 3 Gyr (z = 1.5) to 6 Gyr 
(z — 0.8) was derived using the burst/elliptical models. We 
set the age of the galaxies to be 3 Gyr at z — 1.5 and then 
evolved them in our chosen cosmology until z = 0.8. Third, 
the transformation from AB magnitudes to Vega magnitudes 
was —0.7 mag. All the sources, LBDS & 3CR, were corrected 
for the (l±z) 4 surface brightness dimming, placing them all 
at a common redshift of 0.8 (the 3CR sample has a narrow 
redshift range of ±0.07, and correspondingly small surface 
brightness corrections of < 0.1 mag). We tested the reliabil- 
ity of the corrections by using the two different sets of mod- 
els (Bruzual/Charlot and Fioc/Rocca-Volmerange), and by 
varying the age of the LBDS galaxies at z — 1.5 by ±1 Gyr, 
corresponding to the range of spe ctroscopic ages de rived by 



ourselv es (e.g., Dunlop et al. 199E ) and others (e.g., Bruzual 
& Magi is 199£). The combined uncertainty in K-corrections 
and evolution was ±0.2-0.3 mag. 

We have recalculated the normalization of the 3CR Ko- 
rmendy relation for our cos mology, fixing the (cosm ologv- 
independent) slope at 3.2 ( McLure fc Dunlop 2000 ), and 
plot the fit in Fig. | (solid line). It is seen that 53W069 
and 53W091 both lie on the 3CR Kormendy relation. The 
effective radii of the two LBDS galaxies (r e ~ 3 kpc) are 
much smaller than the mean effective radius of the 3CR 
sources (f e — 12 kpc), but they are comparable with the 
size of the sm allest 3CR galaxy in this sample (3C340 with 
r e = 4.0 kpc). McLure & Dunlop (2000) demonstrated that 
these 3CR galaxies follow a Kormendy relation that is indis- 
tinguishable from that of AGN hosts at z ~ 0.2, assuming 
passive evolution of their stellar populations. Here we can 
extend that result to higher redshifts, concluding that there 
is no evidence for any significant dynamical evolution of the 
elliptical hosts of radio sources between z ~ 1.5 and 0.2. 
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Figure 6. Effective surface brightnesses of 53W069 & 53W091 
compared with the Kormendy relation of a sample of 3CR radio 
galaxies at z ~ 0.8. 53W069 (square) & 53W091 (circle) have 
been corrected for: (i) surface brightness dimming, (ii) bandpass 
shift (observed F110W to F814W), and (iii) passive evolution, 
between z ~ 1.5 and z = 0.8. The 3CR surface brightness data 
(triangles) have been corrected (< 0.1 mag) to a common redshift 
of 0.8. The solid line is t he best-fitting Kormendy relation from 
McLure & Dunlop (2000), renormalized to a redshift of 0.8 and 



our cosmology. 
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Figure 7. Mean surface brightnesses of 53W069 & 53W091 com- 
pared with the Kormendy relation of a sample of cluster ellipticals 
at 2 ~ 0.4. 53W069 (solid square) & 53W091 (solid circle) have 
been corrected for: (i) surface brightness dimming, (ii) bandpass 
shift (observed F110W to B), and (iii) passive evolution, between 
z ~ 1 .5 and z = 0.0. T he cluster ellipticals have been taken 
from Zieglcr ct ah (1999), and have similarly been transformed 
to zero redshift. The clusters are: Abell 370 at z = 0.375 (tri- 
angles); CI 1447±26 at z = 0.389 (diamonds); CI 0939±47 at 
z = 0.407 (plus symbols); CI 0303±17 at z = 0.416 (crosses); and 
CI 0016±16 at z = 0.55 (squares). The solid line is a least squares 
bisector fit to the cluster data. 
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4.2 Comparison with cluster ellipticals at z ~ 0.4 

In Fig. (?], we compare the LBDS radio galaxies with clus- 
ter ellipticals at z ~ 0. 4. This sample cons ists of all the 
early- type galaxies from Ziegler et al. (1999), drawn from 



five clusters observed with WFPC2 — Abell 370 (z = 
0.375), CI 1447+26 (z = 0.389), CI 0939+47 (z = 0.407), 
CI 0303+17 (z = 0.416) and CI 0016+16 (z = 0.55). 



Ziegler et al. (199!:) corrected their data to zero-redshift 



B-band, and compared the differences in average surface 
brightness (AM) between the z ~ 0.4 clusters and the 
Coma cluster (they used two different samples of Coma 
data a nd two different me asures of average surface bright- 
ness). Ziegler et al. (199E) demonstrated that the slope of 
the Kormendy relation for each of these clusters was con- 
sistent with that of Coma, and that the difference in sur- 
face br 



ghtiiesses AM was due Co passive evolution of Che 



stellar populationo, poooibly with 



low level ongoing 



star formation. In Fig. ^ we plot their fully corrected data 
for our cosmology, using the 'ComaSBD/AMf rcc ' surface 
brightness corrections from their table 4 (taking their al- 
ternative AM values does not change our re sults). Using 
a least-squares bisector fit (Isobe et al. 1990) we find the 
best-fitting Kormendy relation for the combined sample to 
be (fi e ) = (2.8 ± 0.2) logr e + (19.4 ± 0.1). 

For each of the two LBDS radio galaxies, we have mea- 
sured the mean surface brightness within r e from the F110W 
data (Table ^) and corrected these (jj, e ) values to zero red- 
shift. Specifically, we applied corrections of: (i) +1.8 mag 
K-correction from observed F110W to restframe B; (ii) 
+ 1.5 mag of passive evolution from 3 Gyr (z — 1.5) to 
13 Gyr (z — 0); (iii) —0.7 mag AB to Vega systems; and 
(iv) —4.0 mag due to cosmological surface brightness dim- 
ming. The combined error in these corrections is ±0.3 mag, 
noting that most of the uncertainty in the models occurs at 
early ages and thus the error is no larger than that in the 
previous section. 53W069 (solid square) and 53W091 (solid 
circle) are overplotted on the cluster Kormendy relation in 
Fig. [71 The LBDS galaxies lie on the Kormendy relation, en- 
tirely consistent with the scatter in the relation for cluster 
ellipticals. 

Our investigations of the Kormendy relation in- 
dicate that 53W069 & 53W091 are passively-evolving, 
dynami cally- relaxed elliptica ls at z ~ 1.5. However, the re- 
sults of Zieg ler et al. (199s| ) show that early-type galaxies 
with weak disk components and low-level star formation, 
do not differ significantly from ellipticals on the Kormendy 
relation. Therefore the results of this section alone do not 
exclude the possibility of some star formation in the LBDS 
galaxies; indeed the disk component identified in 53W091 
(section 3) may be the site of recent low-level star formation 
activity. 



5 SCALING RELATIONS AND MASSES 

We have shown above that the two z ~ 1.5 LBDS radio 
galaxies 53W069 & 53W091 are significantly smaller than 
the average 3CR galaxy at z ~ 0.8 (section 4.1). In this 
section, we consider whether their smaller size could also 
explain their lower radio luminosity. We recall how both size 
and radio power scale with the mass of the galaxy and its 



central black hole, and then estimate of the masses of these 
two radio galaxies. 

The correlation between the central black hole mass 
(Mbh) and t he spheroid mass (A/ sp h) of galaxies is now well- 
established (Magorrian et al. 1998): 



M bh oc M 



sph 



(1) 



For spirals, M sp h is the mass of the bulge (spheroid) com- 
ponent (i.e. excluding the disk), for ellipticals it is the to- 
tal (virial) mass of the galaxy. Further investigations have 
also revealed correlations between Mbh and spheroid lumi- 
nosity (L) and between Mbh a nd stellar velocity dispersion 
(Kormendy & Gebhardt 2001). In the absence of velocity 



dispersion data we adopt the well-known relation between 
galaxy mass and luminosi ty, using the mass-to-light ratio 
for elliptical galaxies from Jorgensen, Franx & Kjaergaard 



(1996 



M sph cc L (2) 

For a de Vaucouleurs r 1 ' 4 law, the integrated galaxy lumi- 
nosity is given by I a I e r^ where the intensity at the ef- 
fective radius (I e ) is trivially related to the effective surface 
brightness /j, e oc —2.5 logic. We adopt a Kormendy relation 
with a slope of 3, fi e oc 3 log r e , this being co nsistent with 
the 3CR results of 



McLure 



relation found in section 4.2 above and the z 



Dunlop (2000[ ), the cluster 
0.2 AGN 



sample of Dunlop et al. (2002). Combining these equations 
thus gives a scaling relation between the effective radius and 
the spheroid mass of elliptical galaxies: 

,1.0 



oc A L 



(3) 



The error in the index is due to errors in the M sp h-L cor- 
relation (e.g. Kormendy & Gebhardt 2001 find Mbh oc L 1 , 



c.f. equations 1 & 2) and the slope of the Kormendy relation 
(±0.2, see values in section 4). The combination of these er- 
rors leads to an error in the power-law index of no more than 
±0.2. 

Whatever the physical mechanism that drives the radio 
output of AGN, and indeed non-active spheroids, a num- 
ber of authors have found a power-law dependence between 
radio power (P) and black hole mass: 

(4) 



P cc M b \ oc M s ; h 



where we have u sed the Magorrian relation (equation 1) in 
the second step. Franceschini, Vercellone & Fabian (199£) 



find a best-fitting value of 7 = 2.7 ± 0.3 for a sample of 
nearby galaxies, whose black hole masses are based directly 
on stellar dynamics. They note that 7 ~ 2 is a generic pre- 
diction of any emission process that depends on the emitting 
area around a black hole, and in particular 7 = 2.2 is the 
prediction of advecti o n-dominated accret ion flows (see e.g., 
Narayan & Yi 1995). Lacy et al. (2001) quote a value of 



7 = 1.4 ± 0.2 for a sample of radio-loud quasars, although 
inspection of their Fig. 2 suggests 7 ~ 2-2.5 is also consis- 
tent with their data. Dunlop et al. (2002 ) demonstrated that 
with 7 = 2.5, equation 4 is consistent with both (i) normal 
galaxies and radio-quiet quasars, and (ii) radio galaxies and 
radio-loud quasars, but with a normalization that differs be- 
tween the two classes of radio power. 

If equations 3 & 4 are combined, we get a scaling rela- 
tion between the radio power and effective radius of radio- 
loud ellipticals: 
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Par 



(5) 



where we have used 7 = 2.7 from Franceschini et al. (1998) 



Given the uncertainties in 7 this can only be considered an 
approximate relation, but it is still interesting to compare 
the two LBDS galaxies with the more powerful 3CR sources 
at z ~ 0.8. We converted the observed 8.4 GHz flux densities 



of the ten 3CR sources (Best, Longair & Rottgering 1997 



to emitted 1.4 GHz luminosites using a mean spectral in- 
dex of ct = 0.8 (where S v tx v~ a )- Their mean power is then 
-P3CR = 9 x 10 27 W Hz -1 in our cosmology. We similarly cal- 
culated the 1.4 GHz luminosites of 53W069 fe 53W091 (with 



spec tral indices of 0.9 fc 1.1 respectively; Waddington et al 



S nnradet al. 1997), giving Plbds ~ 5x10™ W Hz ~" 



Given these radio powers, and recalling that the mean ef- 
fective radius of the 3CR sources is 12 kpc (section 4.1), 
the scaling relation (equation 4) would predict the sizes of 
the LBDS galaxies to be 4 kpc. This is in excellent agree- 
ment with their actual effective radii of 3 kpc. The smaller 
size and lower radio luminosity of these galaxies compared 
with 3CR is consistent with them being less-massive galaxies 
with proportionately less-massive black holes; their smaller 
size cannot be taken as evidence of evolution in the size of 
the largest ellipticals. 

We estimated the virial masses of 53W069 & 53W091 
as follows. Although we do not know the velocity disper- 
sions, we can estimate them from the zero-redshift Funda- 



mental Plane relation (lorgensen et al. 1996). Using the ef- 
fective radius and present-day B-band surface brightness of 
each galaxy derived in section 4 ({fJ,) e = 20.4 for 53W069 
& 20.0 for 53W091), we predict a velocity dispersion of 
a ~ 150 ± 45 km s" 1 for the LBDS galaxies. The present- 
day absolute B-band magnitudes (Vega system) were simi- 
larly calculated from the observed F110W magnitudes and 
the corrections of section 4. We find Mb = —20.2 for 
53W069 and —20.3 for 53W091, correspond ing to luminosi- 



ties of L b = 0.9L* B and l.OL^ respectively ( Madgwick et al i 
200 2|). j jsing the Faber- Jackson rel ation recently measure d 
for Coma and three z ~ 0.4 clusters (Ziegler & Bender 1997) 



we estimate a velocity dispersion of a ~ 180 km s~ for 
the LBDS galaxies, consistent with the more accurate de- 
termination from the zero-redshift Fundamental Plane. The 
virial mass, M v i r ~ 5a 2 r e /G, is inferred to be M V - 1T ~ (0.4- 
1.2) x 10 11 Mq for each of the two galaxies. From equation 
(1), the mass of their central black holes is inferred to be 
(1-3) x 10 8 Mf?), taking the constant of proportionalit y to 
be 0.0025 (fVlcLure & Dunlop 200 1|; iDunlop et al. 2002j). 



6 CONCLUSIONS 

We have used HST observations of the two millijansky ra- 
dio galaxies LBDS 53W069 and 53W091 to investigate their 
optical morphologies at z ~ 1.5. In both F814W (restframe 
(7-band) and F110W (restframe B-band) 53W069 is best de- 
scribed by an elliptical (de Vaucouleurs) model, of effective 
radius C//3 or 3 kpc. The U — B colour gradient indicates 
that the galaxy formed in a rapid burst of star formation at 
high redshift (z > 5) and has a high metallicity (3Zq) stel- 
lar population. In F110W, 53W091 is similarly modelled as 
an elliptical galaxy of effective radius 0'.'3 (3 kpc). At shorter 
wavelengths (F814W), 53W091 is more extended, consistent 
with there being a faint blue disk contributing ~ 20% of the 



flux within r e . The U — B colour of the galaxy within a ra- 
dius of 0'.'4 indicates a high metallicity (3Zq), high-redshift 
(z > 5) burst of star formation, whereas the metallicity falls 
to solar at larger radii. 

Assuming passive evolution of their stellar populations, 
these two elliptical galaxies lie on the B-band Kormendy sur- 
face brightness-effective radius relations of both 3CR radio 
galaxies and cluster ellipticals. Their sizes and radio lumi- 
nosities are consistent with scaling relations applied to the 
3CR radio galaxies, in which both radio power and effective 
radius scale with galaxy mass. 

Our analyses of the restframe UV spectra of these two 
galaxies demonstrated that they were p assively-evolving el- 
lipticals, with ages > 3 Gvr at z ~ 1.5 ( Dunlgjj^ct^ al. 199(: ; 



Spinrad et al. 1997; Dunlop 1999b 



Nolan et al 



2002| ). 



The 

age of the universe at this redshift is 4.5 Gyr in the flat, 
lambda-dominated cosmology that w e assume here, requir - 



ing a formation redshift of Zf > 5. In Peacock et al. (199£) 
we compared the primordial density fluctuation spectrum re- 
quired by the existence of these galaxies with that inferred 
from Lyman-a absorbers and Lyman-break galaxies. Those 
results indicated a similarly high redshift of gravitational 
collapse, z c ~ 6-8. 

Here we have demonstrated that the morphologies and 
internal colours of these galaxies similarly show them to be 
ellipticals whose last major episode of star formation was 
at very high redshift. These results imply that the last ma- 
jor merger was also some while in the past, since such an 
event would produce significant star formation and an irreg- 
ular morphology, for which we see no strong evidence at the 
epoch at which we see them. Although it may be that the 
faint blue disk component in 53W091 is the remains of the 
galaxy's last merger. In this context, we note that 53W091 
has several companions with similar colours, one of which 
has been spe ctroscopically confirmed to have the sam e red- 
shift of 1.55 (|Spinrad et al. 1997|; |Bunker et al. 2002|). It is 



possible, even likely, that one or more of these companions 
will merge with 53W091 in the future. In contrast, 53W069 
has no such companions. 

We began the paper by asking how ellipticals formed - 
via monolithic collapse or hierarchical mergers? We end by 
concluding that we cannot be certain, even in the particu- 
lar case of these two LBDS radio galaxies. The most likely 
interpretation of the data is that they have not evolved sig- 
nificantly since high redshifts (z ~ 5), save for passive stellar 
evolution. However, their formation could still have been via 
mergers at early times, and at least in the case of 53W091, 
it seems probably that further mergers lie ahead. What is 
clear, is that these radio galaxies are quite ordinary ellipti- 
cals, and if they continue to evolve passively until the present 
day they will become typical L* elliptical galaxies. 
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